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1 Introduction 

Solution of the two-body problem has exclusive meaning for decoding the structure of 
interactions between both macroscopic bodies and elementary particles. The unique way 
to determine it by experiment is a study the scattering of particles on each other at various 
energies, theoretical description of which is a problem of two particles. A motion of cosmic 
objects (such as satellites, planets, stars, etc.) or a motion of two charged particles also 
reduces to the two-body problem (classical Kepler-Coulomb problem), which a lot of 
works is dedicated to beginning from Newton's formulation of the World Gravity Law, 
Hooke's formulation of elastic deformations and Coulomb's law of electrostatic interaction 
of charges. Laws of Newton, Hooke and Coulomb made possible to solve a bulk of physical 
problems and give rise to a great many discoveries. The up-to-date enunciating of a 
mathematical problem of two bodies interacting via central potential, one can find in the 
book pp. 

The two-object problem, i.e. determination of their trajectories, can be solved in 
principle, if interaction between them is known. On the other hand, Bertrand sets up an 
inverse problem of determining interaction with respect to known trajectories of motion 
of bodies ( [2j). As it is known, according to the Bertrand's theorem only two types of 
central potentials, of Coulomb and harmonic type, give closed circular or elliptic orbits 
(see, e.g., pj). However Bertrand's problem is solved in the assumption that interaction 
between objects depends only on relative distance between them and is central, meanwhile 
as early as in 19 century it was considered that interaction should depend also from relative 
velocity and acceleration. On this way Weber has developed well-working electrodynamics 
which after creation by Maxwell of the theory of electromagnetic field undeservedly has 
been removed aside. Weber has obtained expression for the force of interaction between 
charged particles depending on the relative velocity and acceleration [I]- [TJ. Special and 
general relativity have given a new push to a consideration of the two-body problem. In 
particular, Sommerfeld has considered the relativistic problem in approximation when 
one of masses is infinitely great [8j. Darwin has obtained expression for interaction force 
between charged particles using retarded potential in non-relativistic approximation of 
relativistic Lagrangian [9]. Later on Darwin's Lagrangian was used in quantum mechanical 
calculations. 

Besides dependence on the relative velocities and accelerations, interaction between 
physical objects depends also on mutual orientation of their spins originally interpreted 
classically as angular velocity of rotation ( [10], p. 123), and now as proper moments of 
momenta. In the bound state the spins of interacting objects are arranged up in definite 
way. For example, it is established that in a deuteron the spins of proton and neutron 
are oriented in one direction. Spins of electrons in Cooper pairs have opposite orientation 
whereas electrons and positrons in positronium can have both the same orientations of 
spins (orthopositronium), and opposite orientations (parapositronium). 

By far not simple question about interaction between objects can be clarified by means 
of studying of the equations of motion taking into account all parameters of objects. To 
clarify all aspects of a problem we will start from the classical equations of motion of a 
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material point with internal degrees of freedom, obtained in [TT]- [13]. The aim of this 
paper is to obtain a system of equations of motion, describing non-relativistic motion of 
two interacting spinning objects. Having non-relativistic classical theory it is easy later 
on to pass on to relativistic and quantum variants, using a correspondence principle. 
Equations, which describe any system of N interacting objects with internal degrees of 
freedom, are considered in §2. These equations are applied to system of two objects in §3. 
Unlike Newton's mechanics here it is supposed that additivity of a momentum does not 
take place, i.e. momentum of the system differs from the sum of momenta of subsystems 
by "the momentum excess" which disappears if interaction does not depend on relative 
speed. Center-of-mass variables and relative variables are introduced in §4, as well as 
relations are established between quantities describing a motion of the system as a whole 
and relative motion of subsystems. The equations of motion of the center of mass of 
the system and of relative motion of subsystems are obtained, which contain undefined 
scalar and pseudo-vector functions describing interaction. For a concrete definition of 
these functions and establishing of the spin equations of motion in §5 the moments of 
momentum of the system and its subsystems are considered. Relations between total 
moments of momentum and angular momenta and spins of the system and its subsystems 
are obtained. Equations of motion of spins are considered in §6. Finally, §7 has to do 
with a question of correlation of the energy of system with energies of its subsystems. 



2 System of N objects with internal degrees of 
freedom 

It is shown in Refs. [11]- [13] that the equation of motion of a mass-point with internal 
degrees of freedom whose position is defined by a radius-vector and which interacts 
with external fields, can be presented in the form 

where K — 1,2, N is subscript labeling a mass-point, 

BU 

p K = m K V K = m m V K - — ^ + [S K x W*] , (2) 

OV k 

is dynamical momentum of i^-th point, vtiqk is a naked mass of the point without of 
taking into account an interaction and interior structure, itlk is an effective mass of the 
point, 

BTJ 

Fk = + [Ck x Va ' ] ' (3) 

is a force acting onto K-th point, 

U K = !7x(«,Ri r ,V jr ,W Jfl W K ,..,wjf ) ) = U 0K - ([R K x V K ] ■ C K ) , (4) 
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is a potential function of K-th. point, which generally may be function of coordinates R^, 
velocities = Il K and accelerations W (m ) = d m W K /dt m , m — 0, 1,2, ...,n, of both 
points. 

A mass-point with internal degrees of freedom should be considered as non-inertial 
extended object with internal structure, defined by pseudo- vectors and Ck, which 
also depend on the interaction of the object with external fields. For the system of 
N objects, interacting both with external fields and with each other, functions Uk and 
pseudo- vectors S^, Ck may be represented as sums 

N 

u K = ur + Y, u TK, (5) 

J=l 

N 

- S 0K + S^* = Sqk + Sq^ + S l j K , (6) 

j=i 

TV 

Ci (~\ext /-i I r~\ext I \ v /~\int (n\ 
OK + ^R - ^OK + ^ 0K + ^JK > K I ) 

J=l 

Here functions Uf^ 1 , S^, Cj^ are specified by interaction of objects with external fields, 
whereas U} n K , S l j K , C v j K are specified by interaction of objects J and K with each other. 
Pseudo- vectors 



J K 



Sext cext , \ ^ q 
K — °0if ~r 2-^ 



N 

int 
JK > 

7=1 



N 
J=l 

may depend on the same variables as potential function (4). Internal parameters Sqk 
and C k are specified exclusively by internal structure of objects and do not depend on 
external variables. In Ref. [13] it is shown that for free objects (mass-points with internal 
degrees of freedom) they are expressed in terms of spin, being integral characteristic of 
this structure, as follows 

(10) 

Cok = — ^ok^OK = —^K^0K S K , (H) 

where is a constant with dimensionality of inverse square of velocity, Qqk is a cyclic 
frequency of Zitterbewegung of i^-th object. If <;k = ~ c ~k ■> where ck is some velocity 
the equation of motion for free K-th object is reduced to non-relativistic limit of Frenkel- 
Mathisson-Weyssenhoff equation ( [H]- [16]). If one adopt $k = + c ~k ■> ^ ne corresponding 
equation will describe a particle with opposite direction of spin, i.e. antiparticle. One 
may assume further that all constants ck are equal among themselves and represent 
the velocity of light. However it should be noticed that constant with dimensionality of 
velocity is proportional to product tqk^ok, where tqk is a radius of Zitterbewegung of 
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K-th object. If c K = r 0K Q 0K , the center of mass of free K-th object in its center-of-inertia 
reference frame moves along circle round the direction of motion of its center of inertia 
with velocity ck- It is reasonable to expect the relations (10)-(11), including the cases 
Sox = (ck = oo) and Cox = (Qqk = 0), will be valid also for interacting objects. 
Along with (1) it is necessary to write down the equations of motion for spins 

ds 

—— = [Q NK {t) x s K ] + m K (t) — ffNif(i)[N}f x sk] + mxit) , (12) 

where Q,-^x{t) are angular velocities of precession of spins round directions of vectors 
Nx(t), which can be the same vector for the system of objects, for example, the vector 
of velocity of the center of inertia of the system. rn^(t) has meaning of the moment of 
force acting onto extended object relative to its center of mass. One can try to define 
the structure of m^(t) from analysis of interaction of internal substance of the extended 
object with external fields, which induces possible movement of this substance inside the 
object. 

A geometrical structure of the object, on the one hand, is specified by distribution of 
its internal substance in some spatial volume V which can change due to the interaction, 
and on the other hand, itself determines an interaction of the object with external fields. 
There are only two variants of consideration of extended object: either as a set (discrete 
or continuous) of structureless mass-points, in the same way as it is usually made in the 
mechanics of absolutely rigid or deformable body, or as a set of mass-points with internal 
degrees of freedom, i.e. similar objects. Then by definition expression for spin of K-th 
object looks like 

N K 
i=l 

for discrete set of points, or 

s K = joK^oK = [p K x ^K{pK)]dV K , (14) 
Jv K (t) 

for continuous set of points. Here and p K are radius- vectors of internal points relative 
to the center of mass of i^-th object, 

„(*> = m f >vf > = m f>^ , (15) 

is momentum of i-th point with effective mass m\ K ^ and velocity v\ K ^ in the case of 
discrete set of structureless points, 

*k(Pk) = ^KiPK) = -^-^ , (16) 

is momentum density of elementary mass d/iK, moving with velocity vk(Pk) inside of 
i^-th object, representable as continuous distribution of structureless points. In the case 



5 



when the object with internal degrees of freedom is a set of the points also endowed 
with internal degrees of freedom, expressions for tt\ k ^ and ■n K (p K ) should have the same 
structure as in (2). This variant unjustifiably complicates the description of extended 
objects, but it may be relevant for the description of larger composite objects. 

Differentiation of (13) and (14) with respect to time and comparison of result with 
(12) will give the possibility to determine the structure of m K (t). However at first it is 
necessary to contemplate the two-object problem, and then the problem of many objects 
with internal degrees of freedom. 



3 System of two objects with internal degrees of 
freedom 

In the system of two interacting mass-points, each of which possesses internal degrees 
of freedom, pseudo-vectors Sk and Ck, K — 1,2, due to (6), (7), (10) and (11) are 
represented like 

Si = <TiSi + Sqi* + S^i* , S2 = <?2 S 2 + Sq2* + S^2* , (17) 
Ci = — <TiQ 01 Si + Cqi* + C21 , C2 = — <?2^02 S 2 + Cq2* + C^2* • (18) 

Momenta of these points with respect to (2), (5), (17) look like 

Pit J Text _|_ TTint) 

Pl = moiVl - { l d ^ 21 ] + [(?isi + Sgf + S*?) x WJ = rmVi , (19) 

fit T Text I Tjint\ 

P 2 = m 02 V 2 - 1 \Z 12 1 + [fes 2 + S e f + S»?) x W 2 ] = m 2 V 2 , (20) 

where m 01 , m 1 and m 2, tn 2 are the naked and effective masses of constituents, respec- 
tively. Forces, acting upon the points, with respect to (3), (5), (18) are 

at t Text — I — Tjint\ 

Fi = - ^ + [(-ftOgiSi + Cgf + Cif) x Vx] , (21) 

fttjjext _|_ jjint\ 

F 2 = ~ %r 2 + K-ft^ + Cg* + C*?) x V 2 ] . (22) 

System as a whole, as well as their constituents, is non-inertial object with internal 
degrees of freedom, whose dynamical momentum should have the structure similar to (2), 
i.e. 

dU 

P = m V - — + [(qs + S ext ) x W] = mV , (23) 
o\ 

where m and m are the naked and effective mass of the system, respectively, and potential 
function U and pseudo-vector S ext have to determined exclusively by interaction of the 
system as a whole with external fields. 
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In Newton's classical mechanics, where potential function does not depend on veloc- 
ities, momentum is additive quantity (see, e.g., [IT]). It means that the momentum of 
the system is determined as sum of momenta of its subsystems P = Pi + P2. However 
experimental data of nuclear and elementary particle physics testify to an absence of 
the momentum additivity. For example, momentum of nucleus considered as a system 
of interacting nucleons does not equal to sum of momenta of nucleons. Due to strong 
interaction arises a mass excess. Therefore the composition law for momenta should be 
written down in the form 

P = Pi + P 2 + AP , (24) 

where AP can be named similar to the mass excess as "momentum excess" . In general 
case we shall consider momentum excess to be caused not only by the interaction of 
constituents of system with each other, but also their interaction with external fields. 
Therefore we suppose 

AP = AT> ext + AP ini . (25) 
Momentum (24) should satisfy to the Second Newton's Law 

dP 

— = F = Fi + F 2 , (26) 

the right hand side of which includes the resultant of all forces, acting onto points of the 
system. This resultant force should be of the same structure as expression (3), i.e. 

BTI 

F = - — + [(-tfJgs + C ext ) x V] . (27) 

All quantities entering into expressions (23), (27), characterize system as a whole and 
can be expressed through corresponding partial quantities. Spin of the system s should 
also satisfy to equation of type (12), or 

ds <Jp(t) _ . . , , . 

it = ^ri p x s i + m < () • (28 > 

Thus, on the one hand, the system of two objects with internal degrees of freedom is 
described by system of six equations (1) and six equations (12), which, on the other hand, 
should be equivalent to three equations (26) and three equations (28). The remaining six 
equations should describe internal movements in the system. These equations we shall 
consider in the separate section. 



4 Relations between quantities, describing 
the system and its subsystems 

Quantities m , m, U, q, Qq, s ; S ext , C ext , entering into (23), (27), may be determined, if 
partial quantities moK, rriK, Uk, Sk, &ok, s k, &oki ^oki entering into (1), are known. 
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To find a connection between all these quantities we transform system of equations (1) 
and (12) in standard way by introducing the relative variables 

r = R 2 - Ri , v = V 2 - Vi , w< m > = W { 2 m) - W<f l) = d m w/dt m , (29) 
and the center-of-mass variables 

= m iRi + m 02 R 2 _ m 01 V 1 + m 02 V 2 (m) _ m iWj m) + m 02 W { 2 m) 

where 

™o = m 01 + m 02 . (31) 

We define also the internal variables, namely, radius-vectors, velocities and accelera- 
tions of points relative to the center of mass 

r K = R K -R, VK = V K -V, w<?> = W&° - W< m > . (32) 

From (29)-(32) we obtain 

r x = R x — R = r , r 2 = R 2 — R = r , (33) 

m m 

m iri + m 02 r 2 = m iRi + m 02 R 2 - m R = , (34) 

as well as corresponding relations for velocities and accelerations. If radius- vector R(£) of 
the center of mass and relative radius- vector r(t) are known, then radius- vectors Rk(0 
of points in absolute reference frame can be easily determined from (33). 

For arbitrary functions /(Ri, R 2 , ...), depending on Ri, R 2 and their derivatives 
with respect to time we have 

_dj_ = moidf_ _d£ = m^ 2 df_ + d£_ 

<9Ri m dH dr ' <9R 2 m <9R dr 

and corresponding derivatives with respect to Vi, V 2 , ... . 

Substitution of quantities R^-(i), Vr-(£), W K (t), obtained from (33) into (24)-(26) 
and comparison of the result with (23), (27), give rise to equations, connecting quantities 
relating to the whole system with partial ones 

dU _ 1 d[m 01 (UZ xt + U^) + m 02 (Ul xt + U'^)} 



<9R m <9R 

d{U% xt + U[f - Uf xt - U l 2 f) 
dr 

+ [(^ lSl + ^ 2 s 2 - q Q 2 s + C ext - Cgf - C 2 f - C e Q f - C\f) x V]- 

■ — [(m 2^i^oiSi - m i<? 2 fio 2 S2 - m 02 C e f - m 02 C 2 i* + m iC^f + m iCi 2 ') x v] , (36) 
m 
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dU 1 d[m m (Ut xt + U%*) +m 02 (Ul xt + E/g*)] 

&v ~ ^ &v + 

p)(Tjext — I — TTint _ Tjext _ Tjint\ 

I u \ u 2 3! 12 1 21 : _ ^p e:E * _ ^p in *_|_ 
+[(?s - sisi - ^ 2 s 2 + S ea;t - Sgf - S*f - S^f - Sf 2 ') x W]+ 

+— [(m 02 ftSi - m i^s 2 + m 02 Sof + m 02 S^* - ra iSgf - moiSjf) x w] . (37) 
mo 

For free system which is not undergo to action of external fields it is necessary to put 
U = U{ xt = Ul xt = 0, A~P ext = 0, S ext = S e xt = Sgf = 0, C ext = Cgf = Clf = 0. It 
is reasonably to admit also that Uffi and depend only on relative variables whereas 
U{ xt and C/f** do not depend on them. Since they are scalar functions depending on 
relative variables we have 

Uif = U[f = U int . (38) 
Then (36), (37) reduce to equations 

[(ft^Si + <r 2 ^ 2 s 2 - <fr 2 8 - C%* - Cf 2 *) x V]- 

— -[(m 2<^oiSi - m 01 ^ 2 02 S2 - m 02 C 2 f + moiCjf) x v] = , (39) 

AP mt = [(qs - ftsi - ^ 2 s 2 - S»? - S\f) x W]+ 

+— [(77io2?iSi - m i<r 2 s 2 + m 02 S^' - m iS£f) x w] , (40) 
'"ii 

where one may get rid of dependence on the center-of-mass variables when assuming 

= ?1 si + ^ 2 s 2 + S*f + Sf 2 * , (41) 

= ^ - ng)si + ft^ - ^)s 2 - ng(S*? + Sj?) . (42) 

These relations follow from the Galileo's relativity principle, according to which equations 
(39), (40) should be covariant relative to Galileo's transformations, so that in the center- 
of-mass reference frame (V = 0, W = 0) they take the form 

[(m 02 <^oi s i - m i<^02 s 2 - rn Q2 C^ + m 01 C^) x v] = , (43) 

AP ,ni = — [(m 02ftSl - m 01 ^ 2 s 2 + m 02 S l 2 f - m 01 Sf 2 *) x w] . (44) 

It is naturally to accept equation (43) to be valid in the case of interaction of the system 
with external fields. Then substitution of (41)- (44) in (36), (37) leads to equations 

9U _ 1 d(m 01 UZ xt + m 02 U^ xt ) t t t 

dR~m~ M + i(C " C oi " C 02 ) x Vj + 
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+— [(mo2CS?-moiCS?)xv], 



(45) 



dU 1 d(m 01 U[ xt + m 02 U% xt ) 
dV ~ dV 



+ [(S ext - S e f - S e f) x W] + 



1 



+ — [(moaSgf - moiSg?) x w] - AP e *' . 
mo 

To get rid of relative variables it is sufficient to put 



m 



m 



(46) 



(47) 



where functions or constants S ext , C ext , £ and T are specified by external fields. It follows 
from (47) 



Sext Qext i &ext f~\ext (-text . f~\ext 
— O 01 -t- O 02 , ^ — O 01 -t- ^ 02 

Now equations (45), (46) take simple form 

dU d m m Ul xt + m 02 U% xt 
<9R <9R m 

ext _ ® ( mm~Uf xt + m 02 U% xt 



AP 

It follows from (49) 

l/(t;R,V,...) 



m 



- U 



moi -Ur + ^Ur + u(t;V,W, 



m 



m 



(48) 

(49) 
(50) 

(51) 



where functions C/jf* are equal to U(t; R, V, ...) up to arbitrary function u(t; V, W, ...), 
determined from initial and boundary conditions 



Uf?(t; R, V, ...) = [/(*; R, V, ...) + u(t; V, W, ...) . 
In accordance to (50) it is connected with external momentum excess 

du(t; V, W,...) 



AP ext (t; V,W,...) 



<9V 



(52) 



(53) 



which hence is determined by dependence of initial and boundary conditions from the 
velocity of the center of mass of the system. If such dependence is absent, then we have 
AP ext (t; W,...) = 0. 

Thus, equation of motion (26) of the system as a whole takes the following form 



d_ 

di 



BU 

m V - — + [( ft si + ^ 2 s 2 + S 2 f + S# + S ext ) x W] 



dU 
dR 



- ^[( ftSl + ^ 2 s 2 + Slf + S#) x V] + [C e ^ x V] , 



(54) 
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where functions U = U(t; R, V, ...), S ext = S ext (t; R, V, ...), C ext = C ext (t; R, V, ...), as 
well as S l 12 *(r, v, ...) and S 2 ™'(r, v, ...), specified by both structure of constituents and their 
spins Si and S2, should be determined in advance. 

Let us define now the relative momentum which in view of relations (47), (52) is 



p = P 2 - Px = 



m 02 - mpi 
m 



m V - + [S- x W] 



[fes 2 - ftsi + S™ - S™ ) x W] - 2 



dU 



int 



<9v 



+ 



<9V 

2m m m 02 



+ 



rrin 



(m v+ [S ext x w])+ (55) 



+<;[s x w] [(moiCiSi + m 02 ^ 2 s 2 + m iS 2 f + m iS l 12 ') x w] 

'"n 



and due to (1) and (19)-(22) satisfies to equation 

F 2 — Fi 



dp 

~dt 



(56) 



Thus, system of six equations (54) and (56) describes six translational degrees of 
freedom of the system. 

For free system (U± xt = U% xt = 0, C ext = 0, S ext = 0) the relative momentum equals 

P = (m 02 - m 01 )V + [fes 2 - flSl + Sj? - S«?) x W] + 2m ° im ° 2 v _ 



m 

QJjint 1 

-2-^ — + — [(to 02 <TiSi + m ^ 2 s 2 + m 02 S 2 f + m iS^) x w] , 
av m 



(57) 



whence it is obvious that it depends not only on the state of movement of system con- 
stituents, but also on the state of motion of the center of mass of the system. Equations 
of motion (54) and (56) in this case look like 



d_ 

di 



(m V + q[s x W]) + tf}g[s x V] = , 



(58) 



d 
di 

d 
"di 



(m 2 - moi)V + [fes 2 - ?1 si - S|f + S™') x W] + 



2m 01 m 02 ^ 
m 



QTJ'int 1 

2— ?[s x w] + — [(moifasi + S 2 f ) + m 02 fes 2 + S#)) x w] 

av m 



= -2- 



dr 



+ [(ft^ lSl - &%s 2 - C™ + C™) x V] - &*[s x v]+ 



+— [(moi^i^si - C 2 f ) + m 02 (^ 2 ^ 2 s 2 - C™)) x v] . 
m 



(59) 



In (55)-(59) scalar function U mt (r, v, ...) and pseudo-vector functions S^f(r, v, ...), 
S 2 "*(r, v, ...), C5_2*(r, v, ...) and C 2 T f (r, v, ...), satisfying to relations (41), (42), remain 
indeterminate. They may be obtained from additional equations, following from spin 
equations of motion, which will be considered below. 
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5 Moment of momentum of the system 

As it was told above, each of two constituents of the system may be considered either as 
a system of structureless mass-points or as a system of mass-points with internal degrees 
of freedom. In this paragraph we deal with first variant, when for every i-th point in (1) 
it is necessary to put Sj = 0, Q = (i = 1, 2, A^, Ni + 1, N, N = Ni + N 2 is amount 
of points in the whole system, is amount of points in K-th subsystem). Then, if 
potential function Ui depends on the velocity of the point, the momentum (2) and force 
(3) take standard form Pj = mjVj, = —dUi/dHi, where rrii = m\ K \ R« = R^ and 
Vj = V-^ = dH^/dt are effective mass, radius- vector and velocity of i-th mass-point 
of K-th subsystem relative to the origin of coordinates, respectively. 

Both system and its subsystems are characterized by the moments of momentum 
relative to origin J, J 1; J 2 . Then, assuming the moment of momentum to be additive 
quantity, we have 

N N 

J = ^[R, x = m i[ R i x V *l = Ji + J 2 , (60) 
i=\ i=i 

where 

N K 

J K = ^ m f»[R,«xV«]. (61) 

i=l 

In each of two subsystems one may determine the center of mass defined by radius- 
vector 

N K N k 

Rk = —J2 m f^ K) , rn K = Y. ^ > ( 62 ) 
mK i=i i=i 

whereas radius-vector of the center of inertia of the whole system is 

T3 1 T3 1 (r WnW , (2)t>(2)l miRi + m 2 R 2 , 
R = — > rriiRi = — > m ; R + > m ; R = , 63 

m ^-^ m \ ^-^ j m 

i=i \t=i i=i / 

where m = m 1 + m 2 . Differentiation of (62), (63) with respect to time gives corresponding 
relations for velocities and accelerations. 

Let's note here that if potential function U = U(t; R, V, ...) explicitly depends on 
time, then effective masses can also depend on time. Hence it follows from (62), (63) 

, N K , N K 

V«- = — yj m«V« + — JJ m r»Rf - — — H-k , (64) 
i=i i=i 

v _ miVi + m 2 V 2 + rhiRi + m 2 R 2 _ m R 
m m m 

If we introduce the relative coordinates of the center of inertia of second subsystem 
relative to the center of inertia of first subsystem 

r 21 = R 2 - R x , (66) 
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then (63) gives 

R 1 = R-^r 21 , R 2 = R+^r 21 . (67) 
m m 

Substituting (67) in (65) we obtain 

miVi + m 2 V 2 m x m 2 ~m x m 2 , coN 
V = + r 21 . (68) 

Let be a radius-vector of i-th point of K-Xh. subsystem relative to its center of 
mass Mx, be a radius- vector of the same point relative to the center of mass M of the 
whole system, tk be a radius-vector of the center of mass of fT-th subsystem relative 
to the center of mass M of the whole system, which is determined similar to (62) 



j N K 

r K = m i K)r i > (69) 



Then we have following geometric relations 

Rf } =K + v l = R K + r\ K) , (70) 

r t = r K + r<*> , (71) 

R^ = R + r^. (72) 

Differentiation of (69)-(72) with respect to time gives corresponding relations for velocities 
and accelerations. 

Substituting (70) into (63) we obtain 

N Ni N 2 

^ rriiYi = ^ mf ] Yi + ^ mf ] Yi = miV 1 + m 2 r 2 = , (73) 

i=l i=l i=l 

whence it follows 

mim 2 - mim 2 

miVi + m 2 v 2 H r 2 i = . (74) 

m 2 

Substitution of (71) into (69) gives respectively 

N K N k 

£ "W = o • EK^f* + *. <K,r . m ) = ■ < 75 > 

i=l i=l 

Now taking into account relations (70) and (75) expressions (61) for partial moments 
of momentum look like 

3k = Jmk + s k = L k + J£* + s K , (76) 

where 

n k 

3 MK = L K + J™ = m K [R K x V K ] + [R K xJ2 m[ K) v\ K) ] , (77) 

i=i 
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is a moment of momentum of the center of mass of K-th subsystem relative to the origin, 

N K 

s* = £mWxv(*>], (78) 
i=i 

is proper moment of momentum (spin) of K-th subsystem (i.e. total moment of momen- 
tum of all points of K-th subsystem relative to its center of mass, see formula (13)), 

N K 

J'-^t^x^mf'vf'], (79) 
i=i 

is a moment of internal momentum of K-th subsystem relative to the origin of coordinates, 
which according to (75) vanishes, if effective masses of points of subsystem do not depend 
on time. 

Substitution of (72) into (77) gives 

3 MK = L K + L MK + m K [r K x V] , (80) 

where 

N K 

L K = m K [R x V K ] + [R x £m<* \f } ] , (81) 

i=i 

is orbital angular momentum of K-th. subsystem relative to the origin, 

l mk = m K[r K x v K ] , (82) 

N K 

L MK = L^ K -[v K xJ2rhf ) r ( i K) ], (83) 

i=i 

are orbital angular momenta of K-th subsystem relative to the center of mass M of the 
whole system with and without taking into account of dependence of effective masses on 
time, respectively. 

Substitution of (76) and (80) into (60) and using relations (73), (74), gives rise to 

J = Jmi + Si + Jm2 + s 2 = Li + L 2 + L M i + Lm2 + s x + s 2 . (84) 

Quantity 

/ m n 2 \ 

L = Lx + L 2 = m[R x V] + [R x J] mf'vf' + ^ mfK^ (85) 

\i=l i=l / 

represents orbital angular momentum of the system relative to origin of coordinates O, 
whereas 

N 

s = ^[ri x Vj] = Si + s 2 + L M i + L M2 = Si + s 2 + m 1 [r 1 x Vi] + m 2 [r 2 x v 2 ] (86) 
i=i 
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is spin, or proper moment of momentum of the system, i.e. total moment of momentum 
of all points of the system relative to its center of mass M. 

Relative coordinates and velocity of the center of mass of J-th subsystem relative to 
the center of mass of K-th subsystem are 

r jk = r j - r K = Rj - R K , v JK = vj - v K = Vj - Y K . (87) 

Taking into account relations (73), (74), we obtain from (86) finally 

mim 2 r . , . 

s = si + s 2 H r 21 x v 2 i . (88) 

m 

Comparison of this expression with relation (41) is reduced to relations 

Sfi + S# = (<~ - ft)si + (<" - ^ 2 )s 2 + <~^^[r 21 x v 21 ] , (89) 

m 

CS* + Cj? = (ftflg! - ?ng)si + fe^ 2 - f ng)s2 - ^ im2Q ° [r 21 x v 2 i] . (90) 

m 

Quantity 

i = r^ [r21 xv 2 j (9i) 

m 

is an orbital angular momentum characterizing the relative motion of subsystems. 

Expressions (89)-(91) concern to both objects and are symmetric relative to their 
permutation. Therefore S l j^ and Cjk may be represented as sums of symmetric and 
antisymmetric terms 

Sjt = \<Sjk - Sj)sj + \<sj K - <;k)s k + ^jK^^irjK x v JX ] + S]f K] , (92) 



2 V y 2 V 7 2m 



1 1 2 2 

_ qjK WK*%JK [TjK x Vjk] + (93) 

where mjK = mj + rriK, and are constants associated both with K-th. subsystem 
and system composed from K-th. and J-th subsystem, respectively. 

Let's try to determine antisymmetric terms Sjj^-j and Cj™^ basing on following argu- 
ments. In the process of evolution the state of system changes from some initial state to 
finite one, in which spins of subsystems sk are oriented relative to each other by some def- 
inite way. For example, s 2 = +Si, i.e. As = s 2 — si = 0, in finite state of electron-positron 
system corresponds to orthopositronium, and at s 2 = — s 1 , i.e. As = s 2 — si = — 2si, the 
finite state is parapositronium. Thus, variation of difference Asjx = — As^j = sj — sk 
with time characterizes variation of relative orientation of spin of constituents. From (69), 
(78), (83) we obtain 

Nj N K 
i=l i=l 
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Nj N K 



^mJ J) [ri x Vj] - J^mf ^ x V;] - mj[rj x Vj] + m^[r x x v K ] = 



i=i 



where 



• • T (o) , T (0) 



z=l 



(94) 



(95) 



is a moment of momentum of i^-th subsystem relative to the center of mass M of the 
whole system. 

We find from (88), (91), (94) 



Si 



s - 1 - As s - 1 + As 

o , S 2 = n 



(96) 



Now, writing down the equation (56) in terms of the relative variables and the center- 
of-mass variables and assuming it to be covariant under Galileo transformations, one may 
come to conclusion, that following relations 



d 
~dt 



du 



+ q[s x W] 



dV 

= ^ 2 s 2 - flSl ) = ^As + ^(s - 1) 



+ <&l[s x V] = , 



(97) 
(98) 



M21] 



1 



-(^1^01 S 1 ~ ^2^02 S 2 
0J ^A As f 



^1^01 ~~ ^2^02 ^ g _ ^ 



(99) 



4 4 
should be fulfilled. 

Then the equations of motion (54), (56), describing six translational degrees of free- 
dom, will take the following form 

§ + H + «*[.xVl-[C-xVl = 



d 
~dt 

d_ 

di 



m V - + [S- x W] 



m 01 m 02 

\ 

m 

Qjjint 

+ 



<9V 

Qjjint 



811 

+ l C£Xt x V] = , 



(100) 



+ -q[s X w] + 

ov 4 m„ 



moim ° 2 'S eit xwl 



+ 



"""^-[C^xvl + ^isxv] = 0. 



(101) 



dr m 

where the momentum of the system is given by expression (23). Equation (100) shows, 
that the system in question moves in such a way as if it has no internal degrees of freedom, 
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but interaction with external fields is specified not by potential function U, but function 
U + u. 

Substitution of relations (88), (92), (98) into (55) gives following expression for the 
relative momentum 



m 02 - m 01 

P = 

m 



mo V-^±^ + [S-xW] 



int 



2 dU 
dv 



2 "' ( " " t02 (m v + [S ext x w]) + ^[s x w] = m 2 V C2 - mi V cl . (102) 



m\ ' L " " J/ ' 2 

It follows from here that the relative momentum depends on the state of motion not only 
of constituents but also of its center of mass. 

For a complete solution of the two-body problem the spin equations of motion should 
be added to equations (100), (101) that will be considered in the next section. 

6 Spin equations of motion 

Equations of motion for spins may be written down in accordance with (12) as 

ds K 



dt 

ds 
~dl 
dAs 



= [n v (t) x s K ] + mjf(t) , (103) 

[fi v (f) x s] + m(t) , (104) 

. = [fi v (f) x As] + Am(i) , (105) 
dt 

where Am(£) = m 2 — m\ ^ 0, as the relative direction of spins of interacting subsystems 
can change. 

On the other hand, differentiation of spin (88) with respect to time gives 

ds dsi ds 2 d (m 1 m 2 \ 

— = — + — + — r xv, (106) 

dt dt dt dt\ m 1 V ' v 1 

or 

m = mi + m 2 + T - [fi v x 1] • (107) 
dt 

This relation suggests that pseudo-vectors m, m K should have identical structure. Thus 
we assume that 

X 1 d\ 

m = mi Am = m 2 + -Am = — - + [f2 v x 1] • (108) 

2 2 dt 

Hence, equations of motion (104) may be finally written down in the form 

|(s + l) = [ft v x(s + l)], (109) 
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whence it follows that pseudo-vector s + 1 is constant over module. Thus, system of 
equations (100)-(101) and (105), (109) is sufficient for description of all movements of 
the system of two objects with internal degrees of freedom, where angular velocity of 
precession Ov is constant, if dynamical momentum (23) is conserved. 

It is necessary to notice that the equation of motion of spin of the system in the 
form (109) can take place only in the case when in the system one may distinguish two 
subsystems divided by a relative radius-vector r. If the system is imagined as indivisible 
"atomic object, for which introduction of relative variables has no sense, then in the system 
of two vector equations, (100)-(101), and two pseudo-vector equations, (105), (109), two 
equations, (100) and (109), (with m(£) = 0), are independent, whereas equations (101) 
and (105) lose meaning because of As = 0, U mt = 0, u(t; V, ...) = 0, m 02 = m 01 , v = 0, 
w = 0. In this case spin equation of motion (107) takes the standard form 

f t = [flv(t) x s] . (110) 

The similar situation takes place for the object represented as a set of noninteracting 
mass-points (with or without of internal degrees of freedom). The equations of motion of 
such object are (100) and (110). 



7 Energy of system 

In modern physics the conservation energy law has fundamental meaning. It is shown 
in [TTJ- JT3] that possible explicit dependence of potential function on time and acceler- 
ations of the higher order leads to violation of this law for separately taken non-inertial 
object even if its internal degrees of freedom are not considered. Let us consider a problem 
of the energy of the system in question in detail. 

If the energy conservation law takes place, but is broken for separately taken object, 
it means that the object in question is an open system interacting with its environment. 
Increment of energy of the object is compensated by decrease of energy of external medium 
so that total energy of object and medium remains constant. These general intuitive 
reasoning may be illustrated by example of two interacting objects with internal degrees 
of freedom. 

For one of two objects, both environment and second object are external. Equations 
of motion for objects are given by (1), where expressions for momenta of subsystems 
and forces acting on them, are obtained by substitution of expressions (17)-(18), (47), 
(52)-(53), (92)-(93), (98)-(99) into (19)-(22). It gives 

Pk = moKVK - 9{U + 1 + ^ + [(^s + ^S-) x W K ] = m K Y K , (111) 
Fk- ^ + [( ^7 C -g^os) XV*]. (112) 
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The equation of energy balance, which energy conservation as a special case follows 
from, is due to equations of motion (1) which for the system of two objects with interior 
degrees of freedom are reduced to equations (97), (100) and (101). All these equations 
have identical structure and lead to following equations of energy balance. 

It follows from (1) 

*e k ai u + u + u~) + £ ( 8(u + u+p w; , + 1)) (U3) 



where 



dt dt ^ dW (p 



E K = ^4 + (V, • [is + ^S-) x W K ]- 
2 2 mo 



-(V K ■ V " ) + U + u + U mt , (114) 



d(U + u + U int • 
dW K 

is total energy of K-tYi subsystem in absolute reference frame. On account of relations 
(33), (35) the energy (114) may be represented as 



Er = E MK + E 0K , (115) 

where 

EuK = JSS^l + (v • [(i,s + ^S»<) x W] - (V ■ 8(U + u + U'-) 

m d\ K 



+U + u + m 0K (V ■ v K ) + (v K • [(i ? s + ^S ext ) x W]+ 

2 ttiq 

+(v ' [( 2^ S + ^7 S } x w * ] - ^7 (v - • ~^r~ ) ' (116) 

^ox = ^fi + (vk • [(^s + ^S-) x w £ ] - (-l)*(v* • ^) + E/** , (117) 

is energy of K-th subsystem in the reference frame of the center of mass of the whole 
system. It is obvious that in the center-of-mass reference frame (V = 0, W = 0, ...) the 
energy of K-th subsystem is Ex = E Q k- 

On the other hand, equations (97), (100) and (101) lead to 

dE,, du „ du 



- = -— - V (— — ■ W (k+1) ) (118) 



dt ' dt ^ K dWW 

k=0 

dE dK dU " JU_ {k+l) 

dt dt + dt + 2^y dW ( k) w ) > ^) 

k=0 

^ = ^ + f (^.w(^)) (120) 

k=0 
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where on account of (47)- (48) 



E = ^ + (V.[S-xW]-(V.^±^) + f / + M = 

~ + (V-fcs+S- i )xW]-(V.^) + f/-E v (121) 
is total energy of the system, 

E r = H^l^ + i (v . [ S x w]) + (v ■ [S«* x w]) - (v • ^) + (122) 

2m 4 m„ " <9v 

is total energy of relative movement in the system, 

^ = ? (V.[sxW] + (V~)-« (123) 

is additional energy arising because of internal degrees of freedom. 

It is not difficult to show that total energy (121) is expressed in terms of E\, E 2) E U) 
and E r in the following way 



where 



E = HxE x + fi 2 E 2 +u-£, (124) 



m 2 0l + m 2 02 

£ = 2m 01 m 02 (U + E T )+mlu ^ int _ m (m 02 - m 01 ) _ dU int | 
mgi + mg 2 * m 2 01 + mg 2 <9v 

^ (V ■ [s x W]) - m ^ m ^ (y . gg_t ") ) + (126) 



2(m^+mg 2 ) ' m (mg 1 + m§ 2 ) <9V 

m 0l m 02 (m 02 - m 01 ) 



[mo(V • v) + (v • [S e:rf x W]) + (V • [S ex * x w])] 



™o( m 01 + m 02) 

In the center-of-mass reference frame of the whole system relation (124) reduces to 



E - u(t; 0, 0, ...) = mEo! + fi 2 E 2 - So , (127) 



where 



m 01 m 02 
Z-o — — 2 — ; y 



' (m 02 - ropi) gg/ + m) ' 



+ 

v=o,w=o,... 



m ° M + . (128) 



For free system ([/ = 0, S earf = 0) equation (119) in the center-of-mass reference frame 
reduces to d(E — u)/dt — 0, that gives E — u = /j,iE 01 + /i 2 i?02 — £o = const. It does not 
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follow from equations of motion (1) that energies (117) will be conserved, as the interaction 
potential U' mt , generally speaking, can depend on time. Then, for example, if dE m /dt > 0, 
i.e. energy Eqi increases, function /i2-Eo2 — Sq should decrease, because Eqi cannot in- 
crease infinitely. Since both subsystems are in equivalent positions, energy Eq2 behaves 
similarly. Thus, energies of interacting objects should oscillate with time. It means that 
potential energy U mt of interaction also should be oscillatory function of time. Thereby 
interaction between objects is reduced to permanent energy exchange. Moreover, pres- 
ence of function u(t; V, W, ...), apparently, allows to describe the high-energy interaction 
considered usually from the point of view of relativistic quantum theories. Determination 
of explicit dependence of interaction energy and an ascertainment of a role of function 
u(t; V, W, ...) requires separate careful consideration. 
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